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AN  EVALUATION  OF  COMMERCIAL  DENSIMETERS  FOR  USE  IN  LNG 


J.  D.  Siegwarth,  B.  A.  Younglove  and  J.  F.  LaBrecque 


The  cryogenic  fluids  density  reference  system  has  been  used  to 
evaluate  three  basic  types  of  densimeters:   the  vibrating  element 
type,  the  dielectric  cell  type  and  the  displacement  type.   These 
meters  were  used  to  measure  densities  in  liquid  methane  and  liquid 
methane  mixtures  with  ethane,  propane,  normal  butane,  and  nitrogen. 
Measurements  were  made  over  the  density  range  from  400  to  480  kg/m3 , 
temperatures  from  108  K  to  130  K,  and  pressures  from  1  to  3  bar.   A 
hundred  measurements  were  made  at  various  densities,  temperatures, 
pressures  and  compositions. 

Key  words:  Cryogenic  densimeters;  density  reference  system;  liquid 
methane;  LNG. 


INTRODUCTION 


The  density  reference  system  (DRS)  provides  a  well-characterized  sample 

of  a  cryogenic  fluid  such  as  liquefied  natural  gas  (LNG)  and  related  mixtures 

in  which  densimeters  can  be  placed  for  evaluation.   The  DRS  measures  liquid 

density,  temperature  and  pressure. 

3 
Densities  varying  from  400  to  480  kg/m   have  been  determined  from  the 

apparent  mass  of  a  silicon  single  crystal  weighed  in  the  liquid  by  an  auto- 
matic balance.   Temperature  is  calculated  from  vapor  pressure  thermometers 
filled  with  methane  whose  pressure  is  measured  with  a  quartz-spiral  bourdon 
gauge . 

Mixtures  were  made  from  laboratory  grade  (99.97  to  99.99  mol  %)  methane, 
ethane,  propane,  n-butane  and  nitrogen. 

Three  basic  types  of  densimeters  were  studied:   the  vibrating  element 
type,  the  dielectric  cell  type  and  the  displacement  type. 

2.   THE  DENSITY  REFERENCE  SYSTEM 

The  density  reference  system  is  designed  to  compare  other  densimeters 
to  an  accurate  Archimedes  type  densimeter  in  an  isothermal  liquid.   The 
accuracy  of  this  reference  densimeter  is  discussed  in  a  separate  report  [1]. 

The  liquid  sample,  test  densimeters,  and  other  measuring  devices  are 
contained  in  a  closed  system,  as  shown  in  figure  1.   An  automatic  balance  is 
in  an  upper  gas  space  which  is  an  integral  part  of  the  sample  space.   From 
the  balance  arm  is  suspended  a  silicon  single  crystal,  whose  apparent  mass 
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when  immersed  in  the  liquid  is  directly  related  to  the  liquid  density.   Two 
vapor  pressure  thermometers  measure  the  liquid  temperature  at  the  top  and 
bottom  of  the  sample  [1].   Pressure  of  the  liquid  sample  is  measured  by  a 
quartz-spiral  bourdon  gauge. 

Cryogenic  temperatures  are  achieved  by  flowing  cold  nitrogen  gas  through 
a  coiled  tube  attached  to  the  outside  of  the  sample  holder.   These  coils 
were  used  to  cool  the  sample,  and  electric  heaters,  also  wound  on  the  sample 
holder,  were  used  to  heat  it  when  changing  temperatures  between  data  points. 
Each  data  point  was  taken  with  the  sample  and  container  very  close  to  isother- 
mal conditions.   During  the  period  of  time  data  was  being  recorded,  the  heat 
capacity  of  the  sample  was  used  to  keep  the  temperature  constant.   The  heat 
leak  to  the  sample  was  minimized  by  surrounding  the  sample  cell  with  a 
nitrogen  gas  cooled  shield.   The  shield  temperature  was  maintained  near 
that  of  the  sample  by  controlling  the  flow  of  the  nitrogen  gas  passing 
through  the  shield  cooling  coil.   To  obtain  temperature  equilibrium  rapidly 
in  the  14  liter  vessel,  a  5  cm  diameter  turbine  pump  is  used  to  mix  the 
sample.   The  sample  fluid  is  pumped  through  a  3.5  cm  diameter  intake  line 
from  the  bottom  of  the  sample  vessel  and  is  exhausted  into  the  top  of  the 
vessel  through  a  discharge  line  of  the  same  size.   This  stirrer  is  capable 
of  mixing  the  fluid  very  rapidly  at  full  speed.   The  speed  of  the  pump  has 
to  be  reduced  considerably  during  the  actual  measurements  since  rapid  cir- 
culation of  fluid  sample  affected  the  readings  of  the  densimeters  under 
test.   A  low  stirring  rate,  which  did  not  affect  the  instruments,  was  used 
during  the  measurement  to  maintain  an  isothermal  sample. 

The  density  variation  in  the  sample  holder  was  evaluated  from  the  vapor 
pressure  thermometer  measurements.   The  temperatures  at  the  top  and  bottom 
of  the  sample  holder  were  computed  from  the  vapor  pressure  equation  of 
Goodwin  [2].   From  the  average  of  the  temperatures,  the  density  of  liquid 
methane  was  calculated  according  to  the  work  of  Haynes  and  Hiza  [3].   The 
Haynes-Hiza  (HH)  measurements  agree  with  the  density  correlations  of  Goodwin 
[2]  to  within  the  estimated  uncertainty  of  both.   Since  the  HH  and  the  DRS 
densimeters  are  both  Archimedes  types  [4],  the  DRS  has  been  compared  to  the 
HH  results  directly.   The  density  differences  between  the  top  and  bottom  of 
the  sample  holder  have  been  less  than  0.04%. 

3.   DENSIMETERS 

3.1   Vibrating  Element  Densimeter 

The  basic  principle  relating  density  to  frequency  is  directly  equivalent 
to  the  classical  mass  and  spring  oscillator.   The  vibrating  member  of  the 
densimeter,  a  rectangular  plate  in  one  instance  and  a  vibrating  cylinder  in 
the  other  instance,  has  a  natural  frequency  determined  by  the  stiffness 


constant  and  apparent  mass  of  the  vibrating  element.   The  apparent  mass  of 
the  element  depends  on  the  density  of  the  fluid  surrounding  it  since  in 
vibration  a  portion  of  this  liquid  is  also  forced  to  move.   The  density,  p, 
is  related  to  the  period  of  oscillation,  t,  by 

p  =  A  +  Bt2  .  (1) 

A  and  B  are  constants  determined  empirically  at  cryogenic  temperatures. 

3.2  Dielectric  Cell  Densimeter 

The  Clausius-Mossotti  relation  [5]  is 

p  =  (e  -  l)/(e  +  2)a,  (2) 

where  e  is  dielectric  constant  and  a  is  the  polarizability ,  generally  taken 
to  be  a  constant,  independent  of  density. 

Over  a  limited  range  of  composition,  however,  the  density  can  be  repre- 
sented by  a  linear  function  of  the  temperature  and  dielectric  constant  [5] 
rather  than  equation  (2).   This  linear  relation  can  be  used  along  the  satur- 
ation line  of  the  fluid.   The  densimeter  studied  in  this  work  uses  a  bi- 
linear relation  to  relate  dielectric  constant  and  temperature  to  density. 

3.3  Displacement  or  Archimedes  Densimeter 

The  operation  of  this  densimeter  is  based  on  the  Archimedes  principle; 
that  is,  an  object  immersed  in  a  still  fluid  experiences  a  lifting  force 
equal  to  the  weight  of  the  fluid  it  displaces. 

The  measurement  equation  for  this  device  is 

p  =  p  (M  -  M  )/M  (3) 

o       a 

where  p   is  the  density  of  the  submerged  object  and  M  its  mass.   M   is  its 
apparent  mass.   The  two  unknowns  p   and  M  are  determined  for  this  instrument 
by  calibration  at  two  or  more  well  defined  densities. 

4.   DENSIMETER  PERFORMANCE 

Each  densimeter  was  compared  to  the  DRS .   As  an  independent  check  on 
the  DRS  measurement,  the  density  of  saturated  liquid  methane  was  calculated 
from  the  average  temperature  and  then  compared  with  the  density  measured 
with  the  DRS.   There  were  six  fillings  of  the  system,  five  using  methane  and 
one  with  an  LNG-like  mixture  without  nitrogen.   Nitrogen  was  later  added  to 
two  of  the  methane  fillings  and  the  LNG  mixture.   These  liquids  are  listed 

4 


in  Table  I.   The  symbols  in  this  table  appear  in  all  the  following  figures 
showing  the  pure  methane  and  mixtures  together. 

The  LNG  mixture  contained  approximately  88-1/2%  CH  ,  6-1/2%  C„H  ,  3% 
C  H   and  2%  n-C.H   . 

Table  I. 

Symbol 

1.  Methane  + 

Methane  + 

2.  Same  filling 
Methane  +  1%  nitrogen       0 

3.  Methane  + 

Methane  + 

4. 

Methane  +  2%  nitrogen      A 


Same  filling 


LNG  □ 

5.  Same  filling 
LNG  +  2%  nitrogen  A 

6 .  Methane  + 


Two  types  of  data  are  considered  in  this  section: 

1.  differences,  X.  -  DRS . , 

l       l 

where  X.  is  a  density  determination  by  one  of  the  test  densi- 
meters and  DRS.  is  the  corresponding  DRS  value,  and 

2.  percentage  difference,  100%  (X . -DRS . ) /DRS . . 
Plots  are  of  the  percentage  difference  data  only. 

Four  plots  of  the  percentage  differences  are  provided  for  each  of  the 
densimeter  comparisons:   methane  data  in  chronological  order  and  versus 
density  (different  symbols  used  for  each  filling) ,  and  all  data  in  chrono- 
logical order  and  versus  density  (symbols  given  in  Table  1) .   The  detailed 
analysis  of  the  data  shown  in  the  figures  is  described  below. 

If  the  DRS  and  the  densimeter  under  test  give  the  same  density  except 

for  random  error,  then  the  differences  X.  -  DRS.  for  repeated  density  meas- 

l       l       c  J 

urements  will  tend  to  cluster  about  zero.   In  this  case  the  densimeter  is 
said  to  be  relatively  accurate  with  respect  to  the  DRS.   However,  experience 
has  shown  that  the  differences  often  show  an  average  value  that  is  not  zero. 
In  addition,  this  average  value  may  be  density  or  temperature  dependent. 
This  average  value  may  also  vary  randomly  with  time  or  from  filling  to 
filling. 

In  order  to  include  all  the  above  possibilities,  the  model  used  for  com- 
paring densimeters  to  the  DRS  is 


Xijk  -  DRSijk   =   m  +  s(Di  -  D)  +  63  +  eijk'  (4) 

where  m  is  the  average  value  of  the  difference  over  the  range  of  densities* 
D.  ,  s  is  the  slope  over  this  range,  D  is  the  average  of  the  D.,  6.  is  a  dif- 

1  !  '    -j 

ferential  shift  occurring  on  the  jth  filling,  and  e. .   is  the  random  error 

1  jk 

occurring  for  the  kth  measurement  of  D.  in  the  jth  filling.   The  6.  can  only 

be  evaluated  from  the  methane  data  since  there  were  no  repeated  fillings  of 
the  other  test  liquids. 

The  scatter  in  the  random  components  e. .,  is  expressed  by  their  standard 

XJK 

deviation.   This  will  be  called  the  within  or  within  filling  standard  devia- 
tion.  If  the  6 .  are  also  considered  to  be  random,  their  standard  deviation 
is  called  the  between  or  between  fillings  standard  deviation.   All  these 
parameters  may  also  be  functions  of  the  liquid  compositions.   For  a  glossary 
of  the  statistical  terms,  see  Appendix  1. 

When  an  estimated  standard  deviation  is  given,  the  degrees  of  freedom 
(d.f.)  are  also  given.   Degrees  of  freedom  is  a  measure  of  the  amount  of 
information  available  for  estimating  a  standard  deviation.   The  d.f.  is 
usually  the  number  of  measurements  less  the  number  of  other  estimated  param- 
eters.  For  example,  for  methane  there  were  53  measurements.   If  m,  s  and 
the  four  differential  shifts  6.  are  estimated,  there  would  be  47  degrees  of 
freedom  remaining  to  estimate  the  within  standard  deviation. 

The  experimental  data  are  fitted  to  the  model  given  by  equation  (4) 
using  the  method  of  least  squares.   The  test  of  this  model  is  if  it  provides 
a  statistically  better  fit  than  if  some  or  all  of  the  unknown  parameters 
were  to  be  taken  as  zero.   The  tests  for  the  presence  of  the  parameters  are 
made  at  the  5%  level  of  significance,  i.e.,  there  is  a  5%  chance  that  the 
test  will  say  to  include  a  parameter  that  is  not  needed  to  explain  the  data. 
At  a  smaller  level  of  significance  there  is  less  chance  of  including  an 
unnecessary  parameter,  but  there  would  also  be  less  chance  of  finding  a 
useful  parameter. 

When  estimates  of  m,  s,  and  6.  are  given,  they  are  usually  given  in 
terms  of  confidence  intervals,  e.g., 

m  +  t- (std.  dev.  of  m) 

where  m  is  the  estimate  for  m  and  t  is  a  percentage  point  of  the  t  distribu- 
tion [6]  depending  on  degrees  of  freedom  and  the  level  of  confidence  desired. 
The  least  squares  procedure  used  to  fit  the  data  to  equation  (4)  supplies 


*   The  true  value  of  D^  is  unknown;  the  DRS  value  is  used  instead.   The  range 
of  the  Dj_  is  large  compared  to  the  random  error  of  the  DRS;  so  the  error 
introduced  by  this  should  not  be  important. 


both  estimates  of  m  and  its  standard  deviation.   A  95%  confidence  interval 
used  in  estimating  the  value  of  a  parameter  means  that  95%  of  the  time  when 
such  a  confidence  interval  is  constructed   it  contains  the  true  value  of  the 
parameter . 

No  total  uncertainty  has  been  derived  here  for  any  of  the  densimeters 
tested   since  some  of  the  instruments  were  not  provided  with  low-temperature 
calibrations.   A  total  measurement  of  uncertainty  would  include  limits  of 
error  associated  with  calibration  procedures,  the  limits  of  random  error 
(precision)  and  the  limits  on  those  systematic  errors  that  result  from 
changes  in  the  conditions  under  which  the  densimeter  operates.   Each  of  the 
densimeters  is  compared  to  the  DRS  and  the  components  of  the  error  are  pre- 
sented; i.e.,  those  associated  with  6.  and  e.  ..  . 

In  the  description  of  each  of  the  following  comparisons,  a  table  is 
provided  which  gives  a  breakdown  by  composition  of  the  statistical  signif- 
icance of  the  average  value,  m,  the  slope,  s,  and  the  largest  6..   A  95% 
confidence  interval  is  given  for  each  significant  parameter.   The  table  gives 
the  within  standard  deviation  (WSD)  as  well  as  the  between  standard  deviation 
(BSD)  for  the  5  ..   It  is  not  clear  that  the  5  .  are  random  for  any  densimeter 
--  in  some  cases  it  is  clear  they  are  not  --  so  the  between  standard  devia- 
tions should  be  viewed  in  this  context. 

4.1   Comparison  of  Haynes-Hiza  and  Density  Reference  System  Density  Deter- 
minations . 

As  discussed  earlier,  the  Haynes-Hiza  temperature-density  relationship 
for  liquid  methane  [3],  which  agrees  within  experimental  uncertainty  with  the 
correlations  of  Goodwin  [2],  has  been  used  to  confirm  the  DRS  density  results, 
The  percent  deviations  of  the  HH  density  from  the  corresponding  DRS  density 
determination  for  the  five  fillings  are  shown  as  a  function  of  the  chrono- 
logical order  of  the  data  in  figure  2.   This  replication  of  the  methane  data 
over  a  number  of  fillings  provides  a  means  of  checking  the  reproducibility  of 
a  densimeter  from  one  filling  to  the  next.   In  this  instance,  the  means  of 

the  five  groups,  m  +  6    are  not  statistically  different.   Figure  3  shows  the 

J 
percent  deviation  as  a  function  of  density  for  the  DRS.   The  figure  shows 

what  appears  to  be  a  second. degree  effect.   While  it  would  not  be  surprising 

to  find  this  kind  of  systematic  difference  between  the  two  methods,  the  data 

does  not  support  this  at  the  5%  level  of  significance.   However,  an  offset  is 

detectable  at  the  5%  level.   See  table  2.   This  offset  is  well  within  the 

systematic  error  of  the  DRS  [1]. 


Table  2.   Statistical  summary  for  HH-DRS 

Methane 

3 
Significant  at  5%  level? 95%  confidence  interval  (kg/m  ) 

m  yes  .024  +  .018  kg/m3 

s  no 

Largest  6 .  no 

3 

WSD   =   0.065  kg/m 

BSD   =   not  significant 

Here,  as  in  all  the  following  comparisons,  a  standard  deviation,  whether 
within  or  between,  combines  the  imprecision  of  both  methods  of  determining 
density. 

Summary 

The  differences  between  the  densities  determined  by  these  two  methods  do 

show  a  small  systematic  offset,  but  no  filling-to-filling  changes,  and  the 

estimate  of  their  combined  standard  deviations  is  0.065  kg/m  .   For  the 

3 
methane  normal  boiling  point  density  of  42  3.6  3  kg/m  ,  this  standard  deviation 

is  equivalent  to  0.016%. 
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4.2   The  Vibrating  Cylinder  Densimeter 

The  vibrating  cylinder  densimeter  is  shown  schematically  in  figure  4. 
The  vibrating  cylinder  secondary  electronic  output  meter  can  be  adjusted  to 
indicate  a  minimum  density  as  the  zero  indication  and  a  maximum  density  as 
the  full  scale  reading.   It  was  found  that  initial  settings  recommended  by 
the  supplier  produced  an  offset  of  several  percent  in  density. 

Since  the  calibration  supplied  with  this  instrument  was  not  suitable, 
the  instrument  was  recalibrated  using  the  experimental  data.   The  resonant 
frequency  of  the  instrument  was  measured  directly  using  a  frequency  counter 
with  an  0.1  Hz  resolution.   Some  of  the  resonant  frequency  and  the  cor- 
responding DRS  density  data  were  used  to  determine  the  constants  of  equation 
(1)  by  means  of  a  least  squares  fit.   The  result  is, 

p   =   -203.3672  +  5719.149/f2  (5) 

3 
where  f  is  the  frequency  in  kHz  and  p  is  in  kg/m  .   This  calibration  equation 

was  fit  to  all  the  data  except  that  of  the  last  methane  run.   The  calibration 
was  done  for  the  purpose  of  comparing  the  instrument  to  the  DRS.   The  analy- 
sis of  the  data  below  indicates  that  a  better  fit  might  be  obtained  if  a 
linear  term  in  1/f  or  temperature  were  included  in  equation  (5) . 

The  percent  differences  for  the  five  methane  runs  are  given  in  figure  5 
in  the  order  taken.   There  are  some  obvious  differences  in  the  operation  of 
the  meter  from  one  filling  to  the  next--the  upward  shift  during  the  second 
day  of  the  first  run  and  the  shift  down  on  the  last  run--that  are  statis- 
tically significant  when  based  on  the  estimated  within  filling  standard 

3 
deviation  of  0.20  kg/m   (47  d.f.).   The  data  for  the  first  filling  was  made 

over  two  days,  and  there  is  a  shift  in  the  data  between  the  two  days,  thus 
inflating  the  standard  deviation  for  that  filling.   If  one  eliminates  the 
data  for  the  first  filling,  the  estimated  within  standard  deviation  is  0.092 
kg/m3  (27  d.f .) . 

Occasionally  a  measurement  by  this  densimeter  indicated  a  considerably 
smaller  density  than  the  corresponding  DRS  measurement.   At  the  same  time  the 
variation  of  the  frequency  reading,  normally  +  0.1  Hz,  increased  to  about  +  1 
Hz.  This  increased  "noise"  in  the  frequency  reading  was  probably  due  to  the 
presence  of  bubbles  on  the  vibrating  element  and  these  data  were  excluded. 
These  noisy  measurements  occurred  after  the  liquid  sample  had  been  rapidly 
cooled,  thus  there  could  be  some  boiling  near  the  vibrating  element.   No 
controlled  tests  were  done  to  see  whether  any  liquid  level  effects  exist.   A 
too  low  liquid  level  could  degrade  both  the  accuracy  and  precision  of  a 
density  measurement. 

Both  the  methane  and  LNG  mixtures  were  used  to  fit  equation  (5) .   The 
estimates  of  m  for  the  various  mixtures,  given  in  Table  3,  indicate  there 

11 


could  be  a  shift  for  the  higher  nitrogen  content;  but  comparable  shifts  did 

occur  for  the  pure  methane  data  (figure  5). 

Equation  (5)  takes  density,  but  not  temperature,  directly  into  account. 

Because  both  methane  and  LNG   mixtures  were  used,  temperature  and  density  are 

not  directly  related;  the  temperature  ranges  for  the  mixtures  are  about  the 

same  in  each  case,  but  the  density  ranges  are  not.   As  can  be  seen  in  figure 

6  for  methane,  the  percentage  differences  show  evidence  of  what  could  be  a 

temperature  dependence,  since  the  measurements  are  at  saturation.   Four  of 

the  values  of  s  for  the  five  fluids  are  significantly  different  from  zero  at 

the  5%  level  (see  table  3),  and  there  is  a  consistency  in  these  estimates 

with  those  for  LNG  being  larger  than  those  for  methane. 

The  LNG  data,  figures  7  and  8,  showed  more  variability  than  the  methane 

3 
data.   The  combined  standard  deviation  for  the  two  LNG  runs  is  0.51  kg/m   (20 

d.f.),  which  is  significantly  different  from  the  estimated  standard  deviation 

of  0.20  kg/m   (47  d.f.)  for  pure  methane. 

The  shifts  seen  in  figure  5  are  significant,  about  0.1%  between  the 

first  and  second  days  of  the  first  filling,  and  0.05%  between  the  last  two 

methane  fillings.   It  cannot  be  said  whether  these  shifts  are  random  in 

nature  or  are  due  to  some  overlooked  systematic  effect. 

Summary 

The  calibration  using  eq .  (1)  removed  most  of  the  systematic  offset  for 
this  densimeter,  but  some  density  dependence  is  indicated.   Whether  the 
shifts  seen  in  the  methane  data  are  random  or  otherwise  cannot  be  determined 
at  this  time.   The  estimated  within  standard  deviation  increases  signif- 
icantly for  LNG. 
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4.3   Vibrating  Plate  Densimeter 

The  vibrating  plate  densimeter  is  shown  schematically  in  fig.  9.   The 
calibration  supplied  with  this  instrument  was  for  more  dense  fluids  at  room 
temperature.   Since  the  direct  density  reading  of  the  instrument  was  1.6% 
high,  the  resonant  frequency  was  measured  directly  using  a  frequency  counter 
with  a  0.1  Hz  resolution.   A  calibration  equation  for  the  vibrating  plate 
densimeter  was  then  determined  by  a  least  squares  fit  of  equation  (1)  to  the 
resonant  frequency  of  the  densimeter  and  the  corresponding  DRS  density.   The 
fit  was  made  to  the  pure  methane  and  the  LNG  mixture  data  only.   No  mixtures 
containing  N-  were  included  in  the  calibration  for  the  reasons  discussed 
below.   The  calibration  equation  derived  is: 

p  =  -  550.1099  +  5430.736/f2,  (6) 

3 

where  f  is  in  kHz  and  p  is  in  kg/m  .   This  equation  was  used  to  compare  the 

density  measurements  of  this  instrument  with  the  DRS. 

Figure  10  shows  the  percent  differences  for  the  five  methane  runs  in 
the  chronological  order  taken.   Figure  11  shows  the  same  data  as  a  function 
of  density.   This  instrument  occasionally  read  a  significantly  lower  density 
than  the  DRS  while  the  random  variation  of  frequency  readings  increased 
significantly.   The  increase  in  noise  was  attributed  to  bubble  formations  on 
the  vane  as  a  result  of  rapid  cooling  of  the  test  liquid  and  the  data  are 
not  shown  in  the  figures. 

There  is  some  shifting  of  the  mean  level  within  and  between  fillings. 

If  the  6 .  are  taken  to  be  random,  the  estimate  of  their  standard  deviation 

J    3 
is  0.20  kg/m   (4  d.f.).   The  data  for  the  fifth  filling  showed  the  largest 

shift.   There  is  no  obvious  explanation  for  these  shifts  in  terms  of  sys- 
tematic effects. 

Figure  12  shows  the  percent  deviation  for  all  liquids  in  chronological 
order.   The  within  run  standard  deviations  among  the  liquids  are  not  judged 
to  be  significantly  different.   See  table  4  for  the  within  standard  devia- 
tions for  each  liquid.   The  estimated  within  standard  deviation  for  the 

3 
differences  using  both  methane  and  mixture  data  is  0.16  kg/m   (88  d.f.) . 

The  mean  values  for  all  liquids  containing  nitrogen  were  found  to  be  signif- 
icantly different  from  zero  in  the  positive  direction  (see  table  4).   The 
LNG  plus  nitrogen  offset  is  statistically  smaller  than  those  for  the  methane 
plus  nitrogen  mixtures.   This  could  have  been  influenced  by  whatever  caused 
the  large  downward  shift  in  the  last  methane  data.   In  regard  to  this,  it 
should  be  noted  that  no  filling-to-filling  effects  were  taken  in  to  account 
in  forming  the  confidence  intervals  for  m  in  table  4. 
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Figure  13  shows  the  percent  deviation  for  all  the  sample  fluids  versus 
density.   Only  the  differences  for  the  LNG  with  no  nitrogen  data  shows  a 
statistically  significant  trend  with  density  and/or  temperature.   There  is 
no  obvious  explanation  for  this,  and  since  it  does  not  occur  for  the  LNG 
data  without  nitrogen,  it  may  be  anomalous.   The  estimated  within  standard 
deviation  does  not  seem  to  change  with  composition. 

Summary 

The  calibration  of  this  densimeter  is  nitrogen  dependent,  but  there 
does  not  appear  to  be  any  strong  evidence  of  a  temperature  or  density  de- 
pendence that  has  not  been  removed  by  the  calibration.   There  are  shifts  in 
the  mean  level  from  filling  to  filling,  the  largest  of  which  occurs  on  the 
fifth  filling.   The  estimated  within  standard  deviation  does  not  show  a 
statistically  significant  change  with  composition. 


20 


Temperature  and 
Pressure 
Transition  Section 


Vibrating  Plate 


VIBRATING  PLATE 
DENSIMETER 


Figure    9. 


21 


QJ 


fO 

s- 


q 


S- 
rc 

E 


CM 

CM 

+ 
CD 

95% 
Confidence 
Interval 

+  1 

4- 

O 

O 

i|9A3[   %g    4P 
}ueoLj.Lu6ts 

to         O 
CD        C 
>> 

C3 

—1 

95% 
Confidence 

Interval 

-.096  +  .08 
.018  +  012 

4- 
T3 

to 

o 

iL3A3[   %g    4p 

to       oo 
CD         CD 

CSJ 

CM 

+ 

95% 
Confidence 
Interval 

+  1 
CM 

<f- 

to 

CTi 

o 

iL3A3L   XS    lp 
;upoij.lu6ls 

to      o 

CD        C 
>> 

CM 

s: 
+ 

95% 
Confidence 

Interval 

00 
O 

+  1 
CO 

4- 

-a 
oo 

CM 

o 

i L3A9L   %g    4p 
^ueoij.  lu6ls 

to         O 
>> 

c_> 

.95% 
Confidence 

Interval 

-.043  +  .047 
-.39  +  .11 

4- 
TJ 

OO 

r-. 

CD 

■o 

•=3- 

O 
CM 

CD 

iL3A3[   ;„g   }e 
q.uPDij.  lu6  j.s 

O        O         00 

C        C         CD 

>> 

s_ 

CD 
+-> 
CD 
E 
ro 

ro 

offset    m 
slope     s 
largest  6- 

Within 

Standard 

Deviation 

Between 

Standard 

Deviation 

>i 

QJ 

X! 

P 
CO 

LO 

V 

QJ 

'd' 

QJ 

UJ 

JZ 

<* 

u 

C 

&H 

p 

0 

T3 

to 

■H 

C 

rrj 

-P 

• 

ro 

QJ 

CO 

>. 

e 

■H 

-P 

01 

> 

■H 

QJ 

>i 

0) 

01 

P 

+j 

T3 

C 

3 

■H 

QJ 

■P 

01 

0 

T3 

X 

C 

si 

•H 

QJ 

H 

4-1 

E 

TD 

O 

C 

0) 

• 

01 

QJ 

£ 

CO 

QJ 

tn 

-P 

OS 

Cn 

O 

Q 

C 

P 

14-1 

(0 

-P 

0 

QJ 

si 

P 

-rH 

c 

00 

+J 

P 

1 

c 

QJ 

QJ 

o 

>1 

> 

C 

•H 

Si 

O 

CO 

-P 

Si 

fO 

TJ 

>! 

-P 

•H 

QJ 

rH 

QJ 

> 

P 

a 

E 

QJ 

3 

a 

T3 

01 

ro 

ro 

T3 

C 

0) 

QJ 

>i 

CO 

X 

E 

QJ 

-P 

si 

QJ 

>, 

-p 

C 

QJ 

-p 

rO 

P 

•H 

QJ 

si 

CO 

01 

u 

■P 

C 

c 

QJ 

01 

QJ 

•H 

E 

QJ 

T3 

01 

rH 

QJ 

X! 

QJ 

p 

X 

• 

c0 

X 

c 

4-1 

01 

+J 

4-> 

QJ 

QJ 

U 

P 

P 

QJ 

s 

P 

O 

3 

01 

0 

QJ 

4H 

-P 

QJ 

p 

a 

X 

Si 

m 

"0 

•P 

-P 

c 

E 

E 

+j 

CO 

\ 

c 

01 

si 

CP 

o 

■H 

0) 

-p 

4-1 

^ 

PI 

C 

p 

LO 

QJ 

p 

QJ 

CO 

QJ 

> 

QJ 

Si 

T 

si 

■H 

TJ 

-p 

-p 

(Ji 

C 

CO 

O 

3 

P 

-P 

p 

01 

0 

QJ 

P 

ro 

o 

ip 

D 

QJ 

E 

o 

.-H 

-P 

S 

"ST 

no 

CO 

QJ 

IT 

E 

> 

& 

,* 

QJ 

\ 

■H 

P 

CP 

QJ 

01 

c 

cO 

.y 

£ 

C 

■rH 

Eh 

QJ 

01 

in 

T3 

c 

QJ 

CO 

aj 

Cn 

•<3- 

QJ 

> 

G 

43 

■H 

CO 

O 

+J 

Cn 

p 

4-1 

o 

4- 

CO 

4-J 

CD 

O 

E 

s_ 

at 

CO 

-o 

c 

S- 

<u 

o 

XJ 

, 

O) 

ra 

+J 

o 

03 

•r— 

o 

i — 

en 

lf> 

CL 

o     • 

i —     CO 

CD 

O    cn 

£Z 

c    c: 

'1— 

O    T- 

+J 

S_  r— 

■ 

03 

-c  i — 

I 

S- 

O    V- 

-Q 

4- 

'r— 

QJ 

> 

ai 

0) 

o 

jr 

4-     i- 

*» 

+J 

O    QJ 
4- 

-C 

C   4- 

cr 

+-> 

O   T- 

LU 

3 

■r-  "O 
+-> 

m 

T3 

C  4- 

^ 

3    •!- 

4-    +-> 

3 

=1 
co 

03    QJ 

2 

03 
<L) 

E 

T3 

CO    ■!— 

03 

CO 

h- 

10 
QJ 

+->    O 

~2L 

■r— 

•i-  .0 

+-> 

co    E 

O 

00 

QJ     co 

CL 

<z 
CD 

T3 

-a 

or   cu 

o 

ai 

Q    S- 

CJ 

IZ 
03 

QJ 
Cn4- 

jr 

£Z   4- 

-M 

■1—     ■!— 

QJ 

XJ    Q 

B 

C 
O 

ai 

Q-      • 

JZ 

CO    +-> 

4J 

aj  c: 
t-    QJ 

4- 

s-   E 

o 

O    QJ 

a 

O    S- 

^_^ 

e 

ZJ 

O 

QJ    l/> 

•i— 

-£=     03 

+J 

+->    QJ 

03 

E 

•I— 

E 

> 

O    QJ 

ai 

S-   -C 

a 

C|-   +-> 

(%)    N0I1VIA3Q 


23 


GD 
C\J 


o 

+  CNJ 


CM 


E 
cm 

>- 

LU 
O 


s- 
<u 
+-> 

CD 

E  +J 
■r-    c 

co  aj 

C     S- 

cu   tu 

-o  4- 

<+- 

CD  -r- 
+->   Q 

fa 

a.    • 

en  -t-> 
c  t- 
■r—  to 
+->   c 

ra     CD 

s-  -a 

-O 

•r-  CD 
>   -C 

+-> 
cu 

jz  4- 
+->  o 


-i->   o 
5  -m 


cd   =3 
s_  4- 


00     fO 

ra 

CD     00 


CO     >,   C7) 
CD    +->    C 


CD  ■!- 

-O   4- 


CD 

-a  ■!- 

^ 

E 

o 

o 

CD 

Q-  >, 

V 

CO    4- 

CD  ■!- 
i-    4-> 

4- 

4_    C 

O 

O  CD 
O   T3 

c 

•r— 

o 

CD 

■r— 

-C     CO 

+-> 

+->   i — 

ra 

O 

•  f— 

E  -a 

> 

O    E 

CD 

s-    >, 

c 

4—    co 



to 

1 — 

CD 

CD 

ro 

(%)    N0I1VIA3Q 


24 


i. 

O) 

+-> 

CD  J_ 
E  CD 
•i-  -o 
CO    s_ 

c   o 

O) 

"O    r— 

cu   o 

to  en 
i—  o 
O-i — 

o 
en  c: 


to  <_> 

s-  cu 

JD  4)i — 

•i-  -C   J3 

>  +->  tO 
+-> 

CD  <+_ 

.c  o   c 


en 

LU 


O 

Q_ 


O   TD 
•i-     CD 


3    CJ  4- 


CD 
rtS   "D 


03     00 
CD     03 


>>   fO 

(O   +-> 
+->  t-     oo 

TO    CO  i — 
XJ     C     O 

CD    XI 
>,TD     E 


Ol 

T3 


Dlh- 


o  c 

r—     CL  CD 

-—    oo  E 

r0     CD  CD 

s_  s- 

<+-     S-  =S 

O     O  CO 


QJ 


o 

S_   4- 
H-    O 


(%)  N0I1VIA3Q 


25 


1 1 1 1 

1         1         1 

,  , 

. 

DENSIMETER 
>              1 

c 

D  CDtD    CD 

LU 
f- 
< 
_l 
0_ 

■ 

CD          C3EJ 
c*J     © 

VIBRATING 

cao     CO 

• 

CO<8^CO 

-«*      ♦                 -»■♦- 

«a           H         ♦   *< 

i.  ♦    +           +  * 

-#- 

o      esc    co 

<a  ««j       +■■*■+  + 

-t*-  ■*■  ■*■     ■*■     + 

■ 

CO  CO     co      co 

-    *  -  + 

^  +      \++ 

1 • 1 • » 

♦ 1 » 

♦ 

o 

o 

lf> 

S- 
QJ    i/> 

-t->    r— 
QJ     O 

e  -q 

1/1     >> 

c:    i/> 
O) 

o 

OJ   I— 

CD 

ro 

** 

a.  >, 

+J 
cn-r- 
c    to 
•i-    iz 
+-*   ai 

(ti  "O 

jd   aj 

■r-    _£Z 

o 

>    -M 

CO 

*— n 

QJ  l+- 

■^ 

rO 

-C    o 

F 

c: 

\ 

-c   o 

cr> 

■r-     4-1 

_*: 

2    O 

' 

•a   3 

>- 

S- 

c_> 

3     ra 

^» 

1— 

to 

ro    to 

V 

CO 

E 

-z. 

ro    -t-> 

LU 

4->     ■!- 

ro    to 

Q 

T3     C 
OJ 

>>-o 

+-> 
■<-    t/5 

c  > 

CNJ 

to  o: 

C  Q 

ai 

CNJ 


OJ 


NT 

CJ 


OJ 


CD 


o   c 

Q.T- 

l/> 

ro 

QJ   T3 

s_   a> 

4- 

S_     •!- 

O 

O   <f- 

O     ■!- 

c 

+J 

o 

QJ     C 

(— 

jz    aj 

+-> 

+->    "O 

OJ 


(%)  N0I1VIA3Q 


26 


4.4   Capacitance  or  Dielectric  Type  Densimeter 

This  capacitance  system  uses  a  relation  linear  in  dielectric  constant, 
e,  and  temperature,  T,  to  determine  p, 

p  =  A  +  B(e-l)  +  CT  (7) 

the  values  of  A,  B,  and  C  having  been  selected  from  the  designer's  studies 
of  properties  of  LNG-type  mixtures  of  methane,  ethane,  propane  and  normal 
butane.   The  dielectric  cell  densimeter  is  shown  schematically  in  figure  14. 

After  the  "first  methane  sample  density  measurements,  this  densimeter 
was  readjusted.   The  data  for  this  first  run  was  then  recalculated  in  ac- 
cordance with  the  adjustment.   The  plot  of  the  percent  deviations  for  all 
the  methane  data  versus  order  --   figure  15  --  shows  significant  shifting  of 
the  mean  level  from  one  filling  to  the  next.   Figure  16  gives  the  percent 
deviations  of  the  methane  data  versus  density.   The  apparent  slope  is 
exaggerated  in  this  case  by  using  percent,  but  even  for  the  density  dif- 
ferences the  slope  is  significant. 

Figure  17  displays  the  percent  difference  versus  order  for  all  liquids 

and  figure  18  shows  the  percent  difference  versus  density.   The  within 

standard  deviation  for  various  liquids  does  not  differ  significantly.   The 

3 
estimated  within  standard  deviation  based  on  all  data  is  0.078  kg/m   (86 

d.f.).   It  will  be  noted — see  table  5--that  the  average  difference  changes 

with  composition.   The  95%  confidence  intervals  given  in  table  5  for  the 

3 
average  differences  are  based  on  the  between  standard  deviation,  0.32  kg/m  , 

and  this  assumes  the  6.  to  be  random.   This  assumption  may  not  be  supported 

by  future  testing.   The  95%  confidence  intervals  for  s  in  table  5  indicate 
that  slope  also  changes  with  composition. 

Summary 

This  densimeter  was  factory  calibrated;  and  reads  high  for  the  liquids 
without  nitrogen,  and  low  for  those  containing  nitrogen.   Of  lesser  prac- 
tical significance  is  a  dependency  on  temperature  or  density.   The  shifts, 
6.,  are  indeed  significant  when  compared  to  the  estimated  within  standard 
deviation . 

Calculations  by  Giarratano  and  Collier  [7]  have  shown  that  the  density 
of  an  LNG  mixture  can  be  calculated  from  the  dielectric  constant  and  tem- 
perature using  eq.  (7)  with  some  error.   The  size  of  this  error  varies  with 
composition.   The  errors  for  LNG  mixtures  containing  no  nitrogen  but  ranging 
from  pure  methane  to  an  LNG  mixture  similar  to  that  studied  here  are  as  much 
as  0.4%.   The  systematic  offsets  for  the  methane  and  LNG  measurements  in 
this  evaluation  show  a  difference  of  0.5%,  in  good  agreement  with  their 
work.   If  nitrogen  is  included  in  the  mixtures,  the  error  range  in  the 

27 


density  determined  from  the  dielectric  constant  and  temperature  becomes  much 
larger.   If  the  nitrogen  concentration  range  is  from  0  to  1.4%,  the  error 
could  be  as  much  as  1.8%.   In  this  experiment,  the  methane  plus  nitrogen 
mixtures  are  nominally  0,  1  and  2%  nitrogen.   From  the  offsets  (figure  17) , 
it  is  estimated  by  interpolation  that  for  a  nitrogen  concentration  range  of 
0  to  1.4%  the  error  could  be  as  large  as  1.9%,  again  in  good  agreement  with 
Giarratano  and  Collier  [7].   To  realize  a  density  measurement  with  a  di- 
electric cell  with  an  accuracy  approaching  the  precision  of  such  a  cell 
requires  that  the  composition,  especially  the  amount  of  nitrogen,  of  the 
fluid  in  question  be  known. 


28 


£f— at 


Mounting  Bracket 


Capacitor  Plates 


Housing 


DIELECTRIC  CELL  DENSIMETER 


Figure    14, 


29 


i_ 

m 

QJ 

.*: 

*-> 

> — 

QJ 

F 

oo 

CC 

(/) 

Q 

cz 

0) 

1 

CJ 

<_> 

CD 

o 

S- 

c 

o 

10 

4- 

+-> 

•F— 

>. 

<_> 

S- 

rrs 

03 

Q. 

E 

OO 

CO 

i — 

CM 

CO 

CD 

-z. 

QJ 

<J  1 — 

s« 

C     03 

+  1 

+  1 

CM 

fc«    QJ     > 
LO   X5    i- 

CM 

'      ' 

+ 

CTi  •!-     QJ 
4-    +J 

O0 
LO 

in 

o 

4- 

CD 

c   c 

T3 

z 

O  >—i 

1--. 

i 

' — " 

_l             <-> 

1 

^ 

lUe3L^LU6LS 

00 
QJ 

00 

QJ 

CM 

^f 

CO 

O 

CO 

o 

' — - 

QJ 

(J   r— 

4- 

C     03 

+  1 

+  1 

CJ3 

6«    0)     > 

XI 

z 

LO   XJ     S- 

oo 

i 

cn  •!-   Qj 
4-   +-> 

lo 

CM 

>tf- 

sz    c 

CM 

CD 

>. — -- 

O  i—i 

• 

c_> 

CM 

' 

LO 

O 

il^l  %SC1G 

00 
QJ 

00 

QJ 

q.uG0j.Lu6i.s 

>> 

>~, 

■=a- 

r^ 

CO 

o 

CO 

o 

CM 

QJ 

z 

O   i— 
C     03 

+  1 

+  1 

"■    * 

5-5 

<f«    QJ    > 

4- 

CM 

ID   XI     S- 

ro 

O-l   -r-     QJ 

CO 

■sl- 

X> 

+ 

4-    4-> 

o 

o 

LO 

r— 

O   i— i 

LO 

> — 

<_) 

<_) 

1 

1 

oo 

O0 

00 

r~- 

iL3A9L  %g   }B 

QJ 

QJ 

CD 

}ubdij.iu6is 

>■> 

>> 

oo 

LO 

CM 

QJ 

CO 

CD 

Z 

U  i — 
C     03 

CO 

CD 

^-^ 

^9 

^S    QJ    > 

r— 

O0   -i-     QJ 

+  | 

+  1 

4- 

+ 

4-    +-> 

c    c: 

CM 

CO 

X> 

r— 

O   ■-< 

r-^ 

CD 

OO 

(_) 

«_) 

i 

■ 

^ 

^3" 

iL3A9L    %S     IB 

00 
QJ 

CO 
QJ 

CO 
O 

^UeDLJ-LUClS 

>> 

>> 

QJ 

oo 

O  >— 

CM 

00 

C     03 

CD 

CD 

oo 

6^    QJ    > 

^3" 

O 

O 

- — - 

1 — 

LO   XJ     i- 

<_> 

O0   -i-     QJ 

4- 

-- — . 

4-    -!-> 

+  | 

+  1 

+  1 

C     C 

X> 

4- 

O   i—i 

LO 

(_> 

OO 

n> 

X! 

CO 

o 

CM 

<3" 

CD 

o 

LO 

' 

<3- 

>* 

1 

1 

CM 

asA3L   %g   }B 

oo 

QJ 

00 
QJ 

OO 
QJ 

CD 

CM 

oo 

q.uB0 lj.  luSls 

>i 

>> 

>> 

+ 

••—> 

5- 

E 

00 

<o 

C 

c 

cu 

-a   o 

-a   o 

+J 

4-> 

S_     •!- 

c   s_  •■- 

CLI 

+-> 

oo 

C    03    +-> 

QJ    03    -M 

E 

QJ 

QJ 

QJ 

•r-  XI     03 

0)  XJ     03 

(0 

OO 

CL 

CD 

-C     C    •«- 

3    C  ■■- 

5- 

4- 

O 

S_ 

+->    03     > 

+->    03    > 

03 

4- 

i — 

03 

•I-    4->      QJ 

QJ    +J     QJ 

CL. 

O 

00 

i — 

3  O0    CD 

CO   O0   CD 

-Q 
03 


QJ 


>1 

0) 

X! 

S-l 
rO 

LO 

T) 

CU 

^r 

0) 

in 

X! 

t 

U 

c 

Eh 

P 

0 

T3 

in 

•H 

C 

(0 

+J 

• 

ro 

QJ 

rfl 

>1 

E 

■H 

4-1 

CO 

> 

•H 

CU 

>. 

CU 

CO 

u 

4-> 

T3 

c 

D 

•H 

CU 

4J 

w 

0) 

T3 

X 

c 

X! 

■H 

CD 

En 

LH 

E 

nd 

0 

c 

cu 

• 

CO 

CU 

x: 

CO 

CU 

Cr> 

-p 

PS 

Cn 

0 

a 

C 

SH 

LLJ 

ra 

4-1 

0 

0) 
X! 

u 

•H 

c 

CO 

-P 

u 

1 

c 

CU 

0) 

0 

>i 

> 

c 

•H 

XI 

0 

rO 

+J 

XJ 

ro 

T! 

>. 

4-1 

■H 

0) 

^H 

0) 

> 

!q 

a 

E 

CU 

0 

a 

T3 

CO 

ra 

n3 

T3 

c 

0) 

qj 

>i 

rO 

£ 

E 

CU 

+J 

X! 

0) 

>1 

4-1 

C 

<D 

-p 

rfl 

^ 

■H 

CU 

X! 

ro 

CI 

O 

4-1 

C 

C 

Q) 

tn 

CU 

•H 

E 

0) 

T3 

CO 

1—1 

0) 

-Q 

CU 

4-1 

x; 

• 

ra 

x: 

C 

4-> 

co 

V 

+j 

CU 

0) 

O 

U 

S-l 

0) 

E 

M 

0 

3 

to 

O 

CU 

4-J 

4-1 

CD 

lH 

a 

X 

x; 

4-1 

oo 

■H 

+j 

C 

E 

E 

4-> 

rO 

\ 

c 

CO 

X! 

rji 

O 

•H 

0) 
-P 

4-1 

X 

2 

J 

c 

U 

LO 

0) 

Jh 

0) 

OO 

0) 

> 

CU 

X! 

"3- 

x: 

•H 

TS 

4-1 

4-1 

Cn 

c 

rO 

0 

3 

U 

4-1 

M 

C/J 

o 

a) 

V-l  oo 

O 

4-1 

D 

01 

E 

o 

•H 

4-> 

\ 

•*3<   OO 

rcj 

CU 

tn 

E 

> 

E 

X 

CU 

\ 

•H 

u 

Do 

0) 

co 

C 

rO 

M 

x: 

C 

•H 

En 

CU 

CO 

LO 

T3 

c 

CU 

CO 

01 

Do 

«tf 

CU 

> 

C 

X! 

•H 

rfl 

O 

p 

CT 

U 

4-1 

30 


cr 

LU 

go 

Z> 

-z. 


O 


U- 

S- 

d) 

O) 

-l-> 

■o 

(U 

i. 

E 

o 

co 

r— 

'Z. 

<B 

<u 

u 

to 

T3 

•r— 

CD 

rn 

C 

Ol 

O 

•(— 

U 

i — 

r— 

c 

o 

i 

m 

c 

■1— 

-l-> 

o 

4- 

•r— 

S- 

u 

sz 

+-> 

fO 

o 

a 

Q. 

0) 

re 

<u 

S_ 

<_> 

.c: 

cu 

_£Z 

c 

>^ 

-t-> 

o 

t|- 

^ 

-!-> 

+-> 

<_> 

!Z 

-a 

c 

<U 

cu 

3 

-o 

s. 

4- 

't- 

3 

in 

<T> 

is) 

TO 

• — 

CD 

(/> 

o 

E 

TS 

_o 

to    >, 
QJ    +-> 


10 


■i-     C     C 

W     II)     (U 

c  -a   s- 

<L>  O) 

■Dlfll- 

o;  4- 

<u  Q  -i- 
C  Q 

ro    en 

E  c  +-> 
o   c 

ai  o.  cu 
-C  oo  e 
4^   ai  cu 

S_  S- 

4-     S-  =3 

O    O  to 

O  fO 

c         ai 

O  <D  E 
•i-    JZ 

+j  +->  cu 

ro  JC 

■r-  e  +-> 

>  O 

a»   s-  «+- 

O   4~    O 


(%)   N0I1VIA3Q 


31 


IX. 

LU 


LU 

Q 


CJ 

< 
Q. 
< 
CJ 


0  o"^ 


.a  ^^sf 


*&**     co  ^3  * 


<o 


O     .0  ^ 


-o 


CD 
CM 


o 


C  -r- 

Ol  . 

C  +J 

+->  to 

•i-  C 

(_)  Ol 

CL 

na  OJ 


CD 


o 


rO 


E   M 


OJ 


Cn 


>- 


LU 
Q 


CD 
OJ 


o 

OJ 


OJ 


CO 


CD 
CD 


CD 

cn 


cn  £.;:£ 


<D 


£Z    i — 

-a  <+- 

oo 

ci: 
o 


OJ 


ra     CD   <D 

x:   c  <+- 

+J  -i-  4- 

C1J    "CI    -r- 

oi  a  >) 
r  w  it- 
+->  ai  -i- 


4-  s-   c 

O  O    QJ 

C  T- 

o  ai 


to 


o 

Si, 


(%)   N0I1VIA3Q 


32 


o 
o 


o 

l_ 

ai 

u- 

s_ 

a> 

4- 

+-> 

O 

QJ 

E 

S- 

cu 

I/) 

T3 

c 

S_ 

OJ 

o 

o 

"O 

CD 

1 — 

CD 

IO 

<_) 

u 

tz 

•  f— 

JO 

CD 

+-> 

o 

•  r— 

1 — 

u 

o 

to 

c 

Q. 

o 

fO 

s_ 

O 

CD 

1 — 1 

or 

qj 

-C     QJ 

c> 

LlI 

+->  i — 

<.o 

CQ 

4-    to 

O    -!-> 

^ 

3 

C     £Z 

3 

<D 

O    -r- 

^ 

S- 
3 

+->  -o 

U    QJ 

co 

C    T- 

(0 

=3   4- 

H 

e 

4-> 

Z 

tO 

fO     C 
QJ 

^^ 

+-> 

CO   X) 

C  1 

O 

fO 

•a 

ra  •!- 

vr 

Q_ 

>,   QJ 

>,+j 

+-> 

■r-     fO 

•r- 

co 

1/1 

sz   co 

C 

QJ  r- 

ai 

■o   o 

T3 

OO    E 

QJ 

ce    >, 

.c 

Q    co 

+-> 

CD  QJ 

1 — 

c  .c 

1 — 

•r-  I— 

C  J 

to 

o      . 

CM 

4- 

o 

CO     C 

£Z 

QJ     QJ 

o 

S-    E 

•  r- 

s_  ai 

+-> 

o    s_ 

<o 

O     Z3 

•  r— 

CO 

> 

QJ     fO 

QJ 

-C    QJ 

a  +->  e 


(%)    N0I1VIA3Q 


33 


— » • 1 

DENSIMETER 

% 

1 *■ 

— i 1      i 

• 1              i 

CAPACITANCE 

m 

dD 
CD 

<b 

A? 

<&> 

*4 

-*►*- 

m 

** 

m 

•tf 

«•* 

m 

«9 

•       $t*~ 

■■ 

' • 1- 

— . 

• »- 

» » • 

CD 


CD 


CO 


CO 


<^> 

o 

LT> 

o 

s_ 

-Q 

oj 

f> 

OJ 

in 

E 

o 

QJ 

CO 

C 

^ 

T3 

OJ 

>, 

O 

4-> 

C 

•r— 

ID 

to 

+J 

c 

•  f— 

OJ 

(_> 

-o 

«3 

c^> 

CL 

OJ 

rO 

-C 

CO 

*■ — * 

U 

+-> 

Vf 

£ 

-EI 

4- 
o 

\ 

+J 

c 

CD 

.c 

o 

_*: 

+-> 

■i— 

-t-> 

3 

>- 

13 

H 

13 

"3 

■^ 

CO 

to 
Ol 

n3 

*z 

E 

LU 

T3 

>1 

Q 

rt3 

i/l 

•o 

OJ 

>>T3 

+J 

CJ 

■  f— 

co      • 

C\J 

c: 

CC    i—. 
Q 

V 

O) 

OJ 

T3 

Olr— 

OJ 

•i-    rcj 

-d 

-o  +-> 

+-> 

o   c 

. — 

Q.-I- 

( — 

l/l 

13 

OJ   T3 
S-     OJ 

C  i 

4- 

S_     ■!- 

O 

O   l+- 

c  > 

O  "i- 

V 

o 

OJ    c 

•  p— 

-£=     QJ 

-t-> 

+->   X) 

n3 

•r- 

•r- 

E 

> 

O    OJ 

OJ 

s-  s_ 

Q 

U-     fO 

00 

c  > 

i — 

CO 

O) 

N~> 

=3 

en 

OJ 


(%)  N0I1VIA3Q 


34 


4.5   The  Displacement  Densimeter 

The  displacement  or  Archimedes  densimeter  is  shown  schematically  in 
figure  19.   This  instrument  consists  of  a  float  that  is  weighed  by  an  elec- 
tronic balance. 

Figure  20  shows  the  percent  deviation  for  the  methane  data  versus  order 
while  figure  21  shows  the  same  data  as  a  function  of  density.   This  instru- 
ment was  supplied  with  a  factory  calibration  that  agreed  with  the  DRS  to 
within  0.15%  as  shown  in  figure  22. 

The  shift  observed  for  the  mean  level  on  the  second  filling  is  signifi- 
cant. The  mean  levels  for  fillings  1,  3,  and  4  (see  table  I)  are  not  dis- 
tinguishable on  a  statistical  basis.   A  95%  confidence  interval  for  the  mean 

3 
level  based  on  these  three  fillings  is  -0.17  +  0.025  kg/m  .   A  95%  confidence 

3 
interval  for  the  average  value  for  filling  2  is  -0.53  +  0.040  kg/m   --  a 

downshift  of  0.36  kg/m  .   There  is  no  obvious  explanation  for  this  shift. 

The  calibration  was  adjusted  for  the  last  filling  to  bring  the  densimeter 

into  agreement  with  the  DRS.   The  mean  of  the  differences  on  this  last 

3 
filling  is  +.038  +  .046  kg/m  which  is  not  distinguishable  from  zero  at  the 

5%  significance  level.   The  percent  difference  for  all  liquids  versus  order 

are  shown  in  figure  22  and  versus  density  in  23.   There  is  an  apparent  shift 

(see  figure  22  and  table  6)  in  the  average  of  the  differences  when  nitrogen 

is  present  in  the  methane,  but  this  is  not  observed  for  the  LNG  mixture.   It 

is  not  possible  to  determine  from  the  data  whether  this  small  shift  (0.025% 

maximum)  is  due  to  the  DRS  or  this  densimeter.   There  are  no  significant 

indications  that  this  meter  is  temperature,  density  or  composition  dependent. 

3 
The  within  standard  deviation  for  all  liquids  is  0.07  kg/m   (91  d.f.). 

Summary 

The  comparison  in  this  case  is  between  two  Archimedes  type  densimeters, 

and  it  was  not  expected  that  composition  would  be  a  factor.   The  means  for 

methane  runs  1,  3,  and  4  and  the  two  LNG  runs  support  this  expectation. 

However,  the  addition  of  nitrogen  to  the  methane  in  both  cases  produced 

detectable  upward  shifts  in  the  means,  though  these  may  as  well  be  due  to  the 

3 
DRS.   The  small  within  standard  deviation  of  0.07  kg/m   (91  d.f.)  make  such 

shifts  easy  to  detect.   There  was  one  notable  and  unexplained  shift  in  the 

3 
mean  level  of  the  differences  —  from  -.18  to  -.53  kg/m   --  between  the  first 

and  second  fillings  of  the  sample  holder.   The  calibration  adjustment  made 
before  the  last  filling  resulted  in  an  offset  that  was  statistically  in- 
distinguishable from  zero.   Temperature  density  or  composition  dependence 
does  not  seem  to  be  a  problem.   The  instrument  calibration  is  such  that  the 
density  measurements  agree  with  the  DRS  to  within  the  stated  accuracy  of  the 
instrument. 
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5.   CONCLUSIONS 

This  series  of  tests  has  shown  that  from  the  standpoint  of  precision, 
even  with  the  run  to  run  offsets  included,  all  the  instruments  are  suitable 
for  low  temperature  density  measurements  on  a  fluid  of  known  composition. 
There  are,  however,  some  problems  with  calibration  as  discussed  in  section 
4.   These  tests  have  not  established  the  stability  of  the  various  instru- 
ments at  cryogenic  temperatures  over  a  long  period  of  use. 

Some  of  the  run-to-run  offsets  of  the  vibrating  element  densimeters  may 
have  been  the  result  of  variations  of  the  fluid  sample  level  which  could  not 
be  determined  accurately.   If  the  level  were  too  low,  the  densimeter  read- 
ings are  definitely  affected  and  some  offset  may  still  exist  at  higher  fill 
levels  since  the  acoustic  mounting  of  the  vibrating  plate  densimeter,  at 
least,  affects  the  performance.   The  calibrations  supplied  with  these  densi- 
meters were  not  suitable  for  use  at  cryogenic  temperatures. 

The  measurements  on  the  dielectric  cell  densimeter  support  the  cal- 
culations of  Giarratano  and  Collier  [7].   Additional  information,  probably 
the  composition,  must  be  an  input  parameter  to  a  density  measurement  if  the 
measurement  accuracy  of  this  instrument  is  to  approach  the  available  pre- 
cision when  the  fluid  composition  varies. 

The  main  obstacle  to  the  commercial  use  of  some  of  these  densimeters  at 
low  temperature  seems  to  be  the  lack  of  a  suitable  calibration.   Also,  a 
means  of  periodically  checking  a  calibration  in  or  near  the  field  site  is 
desirable.   The  best  method  of  calibration  and  field  testing  would  seem  to 
be  through  the  medium  of  a  transfer  standard  rather  than  through  the  prop- 
erties of  a  low-temperature  fluids  which  are  difficult  to  analyze  in  the 
field  and  can  easily  be  contaminated.    The  DRS  could  serve  to  calibrate  and 
maintain  such  transfer  standards. 
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APPENDIX  I 
Glossary  of  Statistical  Terms 

Confidence  interval:   An  interval  in  which  the  true  value  of  an  unknown 

parameter  will  lie  a  certain  percentage  of  the  time  such  an  interval  is 
constructed.   The  percentage  used  indicates  the  degree  of  confidence. 

Between  standard  deviation:   A  measure  of  the  scatter  of  the  random  fluc- 
tuations  that  occurs  in  the  mean  from  one  filling  to  the  next. 

Degrees  of  Freedom:   The  amount  of  information  available  for  estimating  a 

standard  deviation.   It  is  usually  the  number  of  observations  less  the 
number  of  estimated  parameters. 

Level  of  confidence:   A  level  of  probability  at  which  a  stated  hypothesis 

will  be  rejected  on  the  basis  of  the  data  even  though  the  hypotheses  is 
true . 

Within  standard  deviation:   A  measure  of  the  random  scatter  observed  in  data 
taken  under  those  operating  conditions  occurring  within  a  filling. 
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